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ABSTRACT 

o 

^__l We present here a new observational technique, Phase Closure Nulling (PCN), which has the potential to obtain 

(^ very high contrast detection and spectroscopy of faint companions to bright stars. PCN consists in measuring 

CN closure phases of fully resolved objects with a baseline triplet where one of the baselines crosses a null of the 

^—^ object visibility function. For scenes dominated by the presence of a stellar disk, the correlated flux of the star 

|J^ around nulls is essentially canceled out, and in these regions the signature of fainter, unresolved, scene object(s) 

dominates the imaginary part of the visibility in particular the closure phase. We present here the basics of the 
PCN method, the initial proof-of-concept observation, the envisioned science cases and report about the first 
observing campaign made on VLTI/ AMBER and CHARA/MIRC using this technique. 

^1 Keywords: Observing Technique, Optical Interferometry, Complex Visibility, Closure Phase, Spectroscopic 

^H Binaries, Companions 

Dh 1. INTRODUCTION 

O Phase Closure Nulling (hereafter PCN) was introduced in a paper by Chelli et al. (2009) .'■i' In that reference, 

■^^ the authors did not emphasize the very peculiar series of coincidences that lead to the technique described here, 

^ which in retrospect merit to be told. 

It really all stems from a few minutes of observations taken in February 2008 during a night of an AMBER Task 

^-H Force commissioning run of.-^ The purpose of this commissioning run was to find, and if possible to cure, some 

I> obvious instabilities in the transfer function of AMBER that rendered the absolute calibration of the visibilities 

^O very difficult. We found in particular that the entrance polarizers (one per beam) were the cause of a variable 

fringe beating which biased the visibilities. We thus removed the polarizers to test on a series of calibrators that 

the stability of the instrument was indeed improved. We chose a series of calibrators from the ESO calibrator 

list, located at various positions in the sky and also covering a large range of fiuxes. Among these calibrators, 

the brightest was the K5III star a Puppis. 

^D This night run was successful, and demonstrated that the stability of the instrument was greatly improved 

without the polarizers. Then, the polarizers were put back, since otherwise the instrument would not have been, 
in ESO terms, nominal. Finally, the culprit polarizers were replaced by the consortium in Oct, 2008, leading to 
improved stability, as expected. Thus, the observations done during our commissioning run were the one and 
only ever performed on AMBER without polarizers in the Medium-resolution spectral mode at i? = 1500. 

Of the calibrators used, a Puppis stood out nicely, because the closure phase on this star varied from to 180° 
in a few hours: a Puppis, being a bright calibrator was also a large angular-size star and was completely resolved 
by our ^ 90 m longest baseline, as shown in Fig. [T] However, the shallow slope of the closure phase transition 
from to TT radians as the baseline crosses the first null needed an explanation. Part of the explanation was 
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Figure 1. The uncalibrated closure phase of aPuppis (in radians) as a function of spatial frequency. Each of three groups 
of points corresponds to one five minute observation with VLTI/ AMBER at 7? = 1500 in K. Thanks to the spectral 
bandwidth used, each observation spans a range of spatial frequencies, so the transition from to tt radians, as the baseline 
crosses the first null of the stellar disk Bessel function, is relatively well sampled. The data was binned to larcsec"^ 
and the error bars, smaller than the symbols, are plotted. The overall shape is very different from the Heaviside jump 
expected from a stellar disk, with a smooth transition from to n. Explaining this transition led us to develop the theory 
of Phase Closure Nulling. 

quickly found in the literature: crPuppis was since 1907 known as a spectroscopic binary and should never have 
been included in a calibrator list. Thus, a probable reason for the deviation from a uniform stellar disk was the 
presence of the companion. However, crPuppis being of the SBl type, the companion must be quite faint. The 
hypothesis that faint companions could have a visible imprint on the closure phase near the null led us to the 
development of the phase closure nulling concept. 

2. PHASE CLOSURE NULLING: THEORY 

An important property of stars resolved by long baseline intcrferometry is that the coherence of the light decreases 
with increasing spatial frequencies down to zero before increasing again following the well known behavior of the 
Bessel functions. Michelson and PeaseSl used this property to measure the diameter of Betelgeuse for the first 
time. 



2.1 The uniform disk case 

For a star represented by an uniform disk of radius Ri,, the visibility is proportional to a Bessel function of the 
first order, that is 

V4u)^2'-f^. (1) 

As shown in Fig. [21 14 (w) is zero and the correlated flux of the star is canceled out at spatial frequencies 
multiples of 0.61/i?^. Another property is that the phase of the visibility jumps by 180 degrees at each crossing 
of the Bessel function zeros (hereafter called visibility nulls or simply nulls). 
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Figure 2. Sketch of the visibity function for a single star, showing the origin of the n phase shift when one crosses a 
nuU. Left: complex visibility; middle: visibility amplitude; right: phase. Three points, selected at short, middle and high 
spatial frequencies, are represented as colored lozenges. The three corresponding complex visibilities are pictured at the 
bottom as Fresnel vectors in an Argand diagram showing the complex plane. 



2.2 Perturbation by a close faint companion 

If one adds a faint unresolved companion to the previous star, assuming for simplicity that the 
is parallel to the frequency axis, the visibility becomes 



in) = 



„i 2tt us 



system orientation 



(2) 



where r is the flux ratio and s is the separation between the two components. For small flux ratios r, the 
amplitude of the visibility is only slightly modified by the presence of the companion. The stronger effect occurs 
around visibility nulls of the primary where the visibility perturbation is of the order of r. This effect remains 
weak and is beyond the performances of current interferometers as soon as the flux ratio is smaller than 1%. 

More interesting is the phase of the object visibility. The tangent of the phase of the previous system is given by 



tan(j){u) 



r sin(27r us) 



Vi,{u) + r cos(27rus) 



(3) 



As exemplified in Fig.ls] the phase is the result of two contributions: from the star V^{u) and from the companion 
j,g2z7rus_ Except around visibility nulls of the primary, that is most of the time, the companion produces a phase 
signature in the range ±r/14(M), which remains small. However, around visibility nulls, in the frequency ranges 
for which |14(u)| < r, the phase signature of the companion becomes significant, with an exact value of 2TruQS 
at the frequencies uq of the nulls (w ns/Ri, for the first null), much greater than 180 degrees. It follows that, as 
opposed to the visibility amplitude, even for small flux ratios, there is always a frequency interval around nulls 
within which the phase signature of the companion is larger that any systematic error and is thus measurable. 

2.3 The phase closure properties 

Unfortunately, the absolute phase of an interferogram is a quantity difficult to measure as it requires an absolute 
reference that in general does not exist. But, with 3 or more telescopes, one can use instead the closure phase 
(j)c defined as the phase of the bispectrum on 3 baselines /(mi2,W23), with 



^"("12,^23) =< «("12)«("23)**(W1 



>. 



(4) 
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Figure 3. Sketch of the visibility function for a star plus faint companion, with the same layout as Fig. [2] The bottom 
Argand diagram displays the corresponding Fresnel vectors. The imaginary part due to the faint companion induces a 
phase of the visibility that can become quite large near the visibility null at the origin. 
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Figure 4. Closure phase modulo 27r of a double system, formed by an extended uniform disk and a point source, as a 
function of the maximum reduced frequency. The spectral resolution is 7?. = 1500. The 3 frequencies {^23, Ui2,«i3} (see 
Eq. PI have been chosen in the ratio 1,2,3, with the maximum frequency around the first zero of visibility of the primary 
star. The closure phase is displayed for 3 separations (10 i?*, 100 7?*, and 1000 7?*) and a flux ratio of 0.01. The signature 
of the secondary source is dominant in the region around the minimum visibility of the primary (near reduced frequency 
0.61). The thin line with a 180 degree shift at 0.61 corresponds to the closure phase of the primary alone. 



where U12, M23 and M13 are the 3 frequencies transmitted by the interferometer with M13 = U12 + W23, * denotes 
the complex conjugate, and <> represents an ensemble average. The closure phase is a self calibrated observable 
that unlike the phase does not need an absolute reference. In addition, it coincides to the closure phase of the 
observed object, that is 

0c = 0o(wi2) +0o(m23) -0o(wi3), (5) 

where (f>o is the phase of the object spatial Fourier transform. 

Fig. |4] shows the closure phase around the first visibility null of the primary, from a double system formed by 
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Figure 5. Sketch of the setup used at VLTI for PCN observations. The right panel shows the stations of the VLTI and 
the DO-HO-KO triplet used. The left panel shows the (u, v) plane and a 2-diinension map of the amplitude of the Fourier 
transform of a star of 6.4 mas (o-Puppis). The green, red and blue lines correspond to the baselines HO-KO, DO-HO 
and DO-KO respectively, the latter sweeping above the first null (black circle) of the star due to the combination of the 
supersynthesis effect and the spatial resolution spanned by the observational bandwidth. 

an extended uniform disk and an unresolved companion with a flux ratio of 1% at various separations. Note 
the importance of the closure phase signature from the companion. PCN is all about observing closure phases 
near such nulls, where the correlated flux of the star is vanishing, and fitting in the closure phase shape a stellar 
disk and companion, that is, three basic parameters, i?*, r and s. We note that, as shown in Fig. [5] in various 
observational cases with two short baselines and a long one, the closure phase equals the phase of the long 
baseline since the phases on the two short baselines can be approximated by that of an unresolved stellar disk, 
i.e., zero. 

2.4 Performances 

The performances of PCN have been analyzed in ChcUi et al. (2009) .^1 We summarize here the main properties: 

• (pc being proportional to r/V^,, there is no gain in SNR in the null. One still have to integrate enough to 
get the signal out of photon and detector noises. 

• Fibered interferometers have a field of view of one Airy disk. The field of view combined to the necessity to 
resolve the parent star imposes a maximum recoverable separation of Smax/R* ~ B/D, which translates 
to a maximum of ^ 100 i?,* for VLTI. 

• In the photon noise regime, the error on the flux ratio is given by e(r) « Z/ S\K [S being the Strehl ratio 
and K the total number of photons), which is 3 times worse than direct detection. 



2.5 Proof of concept: PCN on a Puppis 

Although sampling quite sparsely the transition region around the first null, the a Puppis observations depicted 
in Figs. [I] and p^ have been interpreted in terms of PCN by Duvert et al. (2010).^ The author best fit was 
-R* = 3.23 ±0.015 mas, a magnitude difference of 5.3 ±0.2 between the giant and its companion, and a projected 
distance at time of observation of 11 mas, that is, 4 stellar radii. With these numbers, assuming that the 
companion is a main-sequence star, they were able to estimate in two independent ways the masses in the 
system, both estimate giving the same result: a ^Mq K5III primary and a A2V 2.2Mq secondary. 
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Figure 6. CHARA/MIRC observation of the star tDra in July 2009, with several crossings of the null with closure phase 
jumps. The plot groups the phase closures on the 4 triplets in all the 24 channels of MIRC/GRISM (H band), as a 
function of the spatial frequency of the longest baseline in the triplet. From left to right: in blue, the triplet E1-E2-W2; 
in green, the E1-E2-W1; in light blue, E2-W1-W2; and in red, E1-W1-W2. 

3. SCIENCE CASES 

There are a number of science cases that can benefit from the possibility to detect, to very high dynamics, the 
presence of a faint companion near a bright star, i.e., inside the airy disk patch of light that a star produces in 
a conventional telescope: 

• PCN can be complementary of spectroscopic diagnostics of binarity, in the cases when the inclination or 
period of the double system become unfavourable- for an example of limitations of radial velocity detection 
of planetary companions). 

• PCN can perform the spectroscopy of the companion, when the S/N per channel is sufficient, the companion 
spectrum being the stellar spectrum multiplied by the flux ratio r estimated for each spectral channel. 

• For already know spectroscopic binaries of type SBl with a giant component that can be spatially resolved 
by the interferometer, a single PCN detection can give the nature of the companion, and subsequently the 
mass of the giant star. The details of the method used are to be found in Duvert et al. (2010).^ 

• PCN can help in finding Brown dwarves companions, and in general extend binarity statistics towards the 
low mass ratio limit. 

• finally, in favourable cases, one could expect to detect Hot Jupiters and at least measure their orbit and 
inclination. 

A rough estimate for the number of targets feasible with the VLTI are of w 100 stars in K and sa 500 in J. 



4. OBSERVATIONAL CAMPAIGN 2009 

We started an observational campaign in 2009 on two interferometric facilities: AMBER/VLTI^ with 3 ATs, in 
Medium-Resolution K {R^ 1500), for which we got 11 nights, and, MIRC/CHARA,'^ inK.R^ 400, for which 
we got 1 useful observing night (see Fig. [6]). 
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Figure 7. First results of PCN on HIP 107315. See text for details. 



In the rest of the section, we focus on the AMBER/VLTI data. The targets were of two types: SBl stars with 
a giant primary, and, giant stars for which no companion is known, which could serve either as a mean to test 
the sensitivity limit, or eventually, would actually possess a faint companion. 

4.1 Observational procedure 

The observational procedure was very simple. We got an estimate of the star diameter using the SearchCallSI 
web service. Given the conservative uncertainties on the stellar diameters obtained this way, the time of crossing 
of the first null was not very precisely known. So we tried to stay as long as possible on the science star, typically 
4 to 5 hours, bracketing the observations with a 1 h observation of two calibrators, usually smaller, thus fainter 
objects, giving much noisier measurements. Normally for closure phase, a calibration is not really necessary. 
However in the case of AMBER, the closure phase obtained on a star will contain the phase closure of the 
so-called P2VM, the calibration matrix used to estimate the coherent fluxes. It is thus necessary to correct this 
closure phase from the one of a calibrator. 

For the data reduction we used the new amdlib3^ algorithm which provides accurate complex coherent fluxes. 

We did not use the VLTI fringe tracker since it operates at H band, and, as can be seen in Fig. [5) if in X the 
stars are fully resolved at the distance of the longest baselines, in H they are already resolved at the distance of 
the second longest baseline, and the small coherent flux at that point prevent a reliable fringe correction: fringe 
trackers should operate at a longer wavelength than the science observation! 

4.2 Preliminary observations 

We focused firstly on the AMBER/VLTI data since we had previous experience of such data with our a Puppis 
observations. Regrettably, AMBER observations in 2009 suffer from an important closure phase noise, typically 
5° rms. This comes from both an integration time per frame (ss 100 ms) that is quite long with respect to the 
coherence time at K, and from a piston-phase dependency induced by a dichroic in the VLTI beams, optical 
element that has been replaced in November 2009. So, for this campaign, we cannot achieve the ultimate in 
sensitivity possible with AMBER/VLTI, and were obliged to develop new techniques to reduce the data. 

We have at the moment run only a general analysis on most of the single stars in our observation list. This 
analysis consists in a two-pass method. We first find the position of the zero of the real part of the bispectrum, 
to get a first good estimate of the stellar diameter. Then, we use this as an initial value to fit R*,s and r. 
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Figure 8. Improved stability of tlie closure phase of AMBER from 2009 to 2010. Plots of the closure phase as a function 
of time for calibrators of different fluxes and at very different elevation in 2009 (left) and in 2010 (right) in respectively 
7h and 3h spans. The stability of AMBER has improved by a factor 50 from a variance of 4.9° in 2009 to a variance of 
0.1° in 2010. 

Amber_MR_Calibrators-2010-01-05A 



-3 




GAM MON - UTOShOO 
HIP43603 - UT08h30 



Difference: variance: 0.20 deg, mean: 0.25 deg 



2080 2100 2120 2140 2160 2180 

Woveiength (nm) 

Figure 9. Improved stability of the closure phase of AMBER in 2010: plot of the closure phase as a function of wavelength 
for 2 calibrators, of different fluxes and at very different elevation, taken at 4 hours intervals. The variance is 0.2°. 

Figure It] show the two main aspects of the process for one of our objects, HIP 107315. In panel (a), we show 
the (m, v) coverage of the observations with the approximate position of the null. One sees that in this case our 
observations cut the null at a constant angle. The right (b) panel show the result of a the fit of i?* , which gives 
extremely precise values for the star diameter in all observations where the null is crossed (filled circles). From 
these observations we derive i?* — 3.81875 ±0.018 mas. In the lower three panels, we plot, from left to right, the 
shape of the closure phase (c), the real part of the bispcctrum (d) and the imaginary part (e), where a sine-like 
signal of very small amplitude can be seen. If this signal comes indeed from a companion (physically linked 
or field star), the separation is not well constrained yet by the observations, but the flux ratio can be reliably 
estimated to be 3.310"'^. Thus, with the 2009 observations that suffer from a large closure phase error, we can 
achieve already detections of faint companions. 

4.3 Update on the AMBER closure phase stability 

The stability of the closure phase of AMBER has improved dramatically since the replacement of the IRIS 
dichroic in the last months of 2009. We made new closure phase observations in early 2010, and, as shown in 
Fig. Is] the closure phase does not vary from more than 0.1° in a few hours (instead of 5° in our 2009 observations). 
This is also visible in the spectral direction, as seen in Fig. |9J where the fixed phase artifacts are stable and time 
and can be calibrated down to 0.2° at least. 



Although not fully understood, the closure phase dependence on the piston induced by the previous IRIS dichroic 



was certainly related to differential polarization effects. In this respect, it is remarkable that our discovery of 
the PCN effect arose from the sole unpolarized observations made with AMBER. 

5. CONCLUSION AND PERSPECTIVES 

Phase Closure Nulling is an observational technique consisting in measuring closure phases of resolved objects 
with a baseline triplet where one of the baselines crosses a null of the object visibility function. In these regions 
the signature of fainter scene objects are seen in the closure phase, and stay out since the correlated flux of 
the star is essentially nulled at these positions. PCN is limited in field of view and needs quite high spectral 
resolution, but can achieve detection of high-contrast objects. We report a first estimate of contrasts down to 
10""^ obtained in 2009 on AMBER with a bad closure phase rms of « 5°. We think that the detection limit can 
be pushed much farther thanks to the the excellent closure phase stability (^ 0.1°) displayed by AMBER since 
a hardware change in VLTI that took place in November 2009. 
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